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ABSTRACT 


This  thesis  theoretically  presents  the  phenomena  in¬ 
volved  in  the  flow  of  an  incompressible  fluid  over  a  wedge 
with  a  second  incompressible,  lighter,  and  less  viscous 
fluid  blown  through  the  surface  of  the  wedge.  A  method  is 
developed  to  determine  the  inner  fluid  layer  thickness,  the 
wall  shear  stress  and  the  resulting  local  drag  reduction. 
The  results  predict  substantial  drag  reduction. 
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I .  INTRODUCTION 


The  ever  increasing  operating  depths  and  speeds  of  modern 
submarines  and  the  much  higher  resistance  to  the  effects  of 
contact  or  nearby  underwater  explosions  all  require  a  dra¬ 
matic  increase  in  the  speed  of  underwater  weapons.  The 
existing  power  plant  systems  seem  to  have  reached  their 
limits  in  power  density  and  radical  changes  in  body  shapes 
are  limited  by  launching  facilities.  The  answer  might  lie 
in  the  field  of  boundary  layer  control.  Possibilities  are 
injection  of  polymers  into  the  boundary  layer  to  maintain 
laminar  flow,  suction  to  remove  potential  turbulence ,  more 
sophisticated  weapon  shapes  to  reduce  form  drag,  or  intro¬ 
ducing  a  gas  layer  between  the  surface  of  the  body  and  the 
liquid. 

Of  all  the  options,  the  two-phase  boundary  layer  seems 
to  be  the  most  promising  possibility  for  obtaining  an  order 
of  magnitude  increase  in  speed. 

In  this  thesis,  a  method  is  developed  to  calculate  the 
properties  of  the  laminar ,  two-phase  flow  produced  by  in¬ 
jecting  a  light,  incompressible  low-viscosity  fluid  into  the 
boundary  layer  of  the  incompressible  flow  over  a  wedge.  This 
is  just  a  step  on  the  long  way  to  a  high  speed  underwater 
weapon. 
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II.  BACKGROUND 


A.  DISCUSSION  OF  RELATED  PAPERS 

An  extensive  literature  search  was  carried  out  and  the 
relevant  papers  are  discussed  below.  The  reduction  of  drag 
by  use  of  polymers  is  deleted  from  this  survey  because  the 
drag  reduction  expected  by  this  method  is  probably  not  suf¬ 
ficient  to  increase  the  maximum  speed  of  the  underwater 
vehicle  substantially.  Amazingly  enough,  the  most  recent 
paper  in  the  field  of  two-fluid  laminar  boundary  layer  flow 
was  published  back  in  1969. 

1.  J.  Pretsch  [1] 

Pretsch  starts  with  Prandtl ' s  boundary  layer  equa¬ 
tions  for  single-phase  flow  with  suction  or  blowings 


and  boundary  conditions  y  *  0:u  ■  0,  v  «  vwj  y/Re  ®  :  u  *  0 
where  u  is  the  velocity  tangential  to  the  surface,  v  is  the 
velocity  normal  to  the  surface,  x  is  the  arc  length  along 
the  surface,  y  is  the  normal  to  the  surface,  U(x)  is  the  po¬ 
tential  velocity,  v  is  kinematic  viscosity,  vw(x)  is  the 

normal  velocity  or  blowing  velocity  at  the  surface  and  Re  is 

Ux 

the  Reynolds  number,  Re  - 
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-«)  the  above  equations 


For  strong  suction  (vQ  -*• 
simplify  to  a  second  order,  ordinary  differential  equation 
of  the  form 


Sfs  +  35  -  0  where  5  -  -  — 
as1  35  v 

and  boundary  conditions  5  ■  0  :  u  =  0 

C  00  :  u  =  U(x) 

The  solution  of  this  equation  is  the  velocity  profile 


For  strong  blowing  (v  ® )  the  assumption  is  made 
that  the  viscosity  term  can  be  neglected  and  the  above  equa¬ 
tions  again  simplify  to  a  second  order,  ordinary  differential 
equation  but  now  nonlinear. 

For  a  number  of  values  for  similarity  parameters, 
Pretsch  gives  analytical  expressions  for  solutions  of  flat 
plate  and  wedge  flow. 

The  assumption  of  negligible  viscosity  is  too  re¬ 
strictive  for  the  case  under  investigation  in  this  thesis. 

2.  E.  J.  Watson  [2] 

This  extensive  report  consists  of  three  parts,  i.e. 
Theory  of  Similar  Velocity  Distributions,  Flow  with  Uniform 
Suction,  and  Flow  with  Variation  of  Suction  Velocity.  The 
report  is  based  on  the  asymptotic  theory  of  boundary  layer 
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flow.  Although  the  author  is  aware  of  the  existence  of 
Pretsch's  work,  no  comparison  of  results  is  made. 

In  Part  I  the  author  arrives  at  a  number  of  tables 
giving  velocity  distribution  and  coefficients  of  the  series 
expansions  as  functions  of  similarity  related  constants  and 
dimensionless  parameters. 

In  Part  II  the  similarity  principle  is  abandoned  but 
the  suction  velocity  is  uniform.  The  velocity  is  expanded 
into  a  series  from  which  series  for  momentum  and  displace- 
ment  thickness  and  skin  friction  are  derived.  In  th*  sub¬ 
sequent  section  a  number  of  applications  are  investigated 
and  the  results  are  totalized  in  a  number  of  tables. 

In  Part  III  the  asymptotic  theory  is  extended  to 
cover  general  two  dimensional  flow  with  arbitrary  distribu¬ 
tions  of  mainstream  velocity  and  suction  velocity.  The 
results  are  again  displayed  in  tables.  The  author  also  de¬ 
rives  the  results  in  Part  I  as  a  special  case  of  Part  III. 

As  this  report  only  deals  with  suction  it  is  not  of 
direct  interest  to  this  study. 

3.  R.  D.  Cess  and  B.  M.  Sparrow  [3] 

This  paper  deals  with  a  two-phase  flow  in  forced 
convection  over  a  flat  plate.  The  inner  fluid  is  obtained 
by  boiling  the  outer  fluid  at  the  wall  so  that  a  vapor  film 
is  formed.  The  authors  arrive  at  the  following  conservation 
equations : 


where  0 


e "  ♦  PR  f.0  •  -  0  where  0  * 

9 

V*  +  fLfL*  *  *  0 
with  boundary  conditional 
plate  surface:  fg(0)*fg*(0)  ■  0,  0(0)  ■  1 


-  t 

w  08 


interface: 


V#>  ■  (tot?)  W1 


V  (0>  ■  V  <  V' 


v,,#» 


Dfvv 


0(nfi)  -  0 

\ 

free  stream:  fL' '  2  as  nL  *►  • 

where  fg  and  fL  are  the  stream  functions  for  respectively 
inner  and  outer  fluid,  T  is  the  local  temperature ,  Tw  is  the 
temperature  at  the  wall,  Tm  is  the  temperature  in  the  free 
stream,  PR  is  the  Prandtl  number,  p  is  the  mass  density  and 
u  is  the  dynamic  viscosity,  f  is  a  function  of  the  dimen- 

9  PC 

s ionless  parameter  ng  defined  by  ng  ■  J  •  ngi  is  the 

value  for  this  parameter  at  the  interface  of  the  two  fluids. 

f_  is  a  function  of  the  dimensionless  parameter  nT  defined 
L  i  ri —  L 


by  n. 


The  subscripts  L  and  g  refer  to  respective 


ly  outer  and  inner  fluid.  The  primes  denote  differentiation 


with  respect  to  the  relevant  dimensionless  parameter.  This 


problem  is  solved  by  using  an  Emmons** Leigh  [4]  table  and  a 
successive  approximation  procedure. 


4.  K.  Nickel  [5] 

In  this  paper,  the  laminar  flow  of  an  incompressible 
fluid  without  heat  transfer  is  investigated  using  a  method 
of  bounding  the  solutions  applied  to  Prandtl's  equations. 
After  introducing  the  well-known  dimensionless  parameters  n 
and  f(n)#  Prandtl's  partial  differential  equations  reduce  to 
the  ordinary,  third  order  non-linear  differential  equation: 

f"  +  ff"  +  8(1  -  f'2)  -  0 

with  boundary  conditions:  f(0)  *  fQ 

f* (0)  *  0 
f'(-)  -  1 

The  solutions  are  reflected  in  a  number  of  graphs, 
especially  "abbildung  8”  giving  f''(0)  as  a  function  of  8 
and  fQ  will  prove  to  be  very  useful  in  this  thesis. 

5.  R.  D.  Cess  and  E.  M.  Sparrow  [6] 

This  paper  is  an  extension  of  the  paper  discussed  in 
Section  3  above  and  deals  with  the  case  of  a  subcooled  fluid. 
To  extend  the  earlier  results  it  was  necessary  to  consider 
the  conservation  of  energy  within  the  fluid.  The  results  are 
again  reflected  in  a  number  of  figures. 

6.  E.  M.  Sparrow.  V.  K.  Jons son  and  E.  R.  G.  Eckert  [7] 
In  this  paper  the  laminar  two-phase  flow,  obtained  by 

injecting  gas  into  a  liquid  at  the  wall  of  a  flat  plate,  was 


investigated.  The  authors  arrive  at  the  following  ordinary 

differential  equations 

for  the  gast  fg' ' '  +  *g*g' 1  *  0 

for  the  liquids  fL"  '  +  fLf L ' '  -  0 


with  boundary  conditions: 

at  the  wall:  fg(0)  ■  *g0  (blowing  parameter) 

f  * (0)  ■  0  (no  slip  condition) 

g 

at  the  interface:  £  (n-4)  *  fT  (0)  -  0  (no  mass  transfer) 

g  gi  L 


R  fL''(0)  (stress  continuous) 


f^' (ngi)  *  f^' (0)  (velocity  continuous) 


for  nL  ♦  •:  fL' (nL)  ^  2  (velocity  approaches  UJ 


The  authors  so*ve  these  equations  by  applying  a 
series  expansion  and  subsequently  a  successive  approximation 
procedure ,  similar  to  the  one  described  in  the  paper  by  Cess 
and  Sparrow  (Section  3  above) . 

The  results  are  given  in  a  table  and  a  graph.  It  is 
proven  that  evaporation  at  the  interface  only  works 


advantageous. 


Under  assumption  of 

-  shear  stress  in  liquid  is  negligible 

-  interfacial  velocity  and  pressure  distribution  are  those 
imposed  by  the  potential  flow 
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-  the  inertia  of  the  gee  is  negligible 

-  the  gas  flow  is  laminar 

-  body  moves  steadily  through  water  at  zero  incidence 

-  gas  and  water  properties  are  constant;  in  particular  the 
density  of  the  gas  remains  constant 

the  problem  of  a  two-phase  flow  over  different  kinds  of 

bodies  is  solved  by  a  finite  differences  method.  Results 

are  given  in  a  number  of  graphs. 

8.  W.  s.  Bradfield,  R.  0.  Barkdoll  and  J.  T.  Byrne  [9] 

This  paper  discusses  experiments  with  a  "gassing" 

model  to  investigate  possible  drag  reductions.  The  vapor 

layers  were  obtained  by  film  boiling,  sublimation  of  a  solid 

surface  and  chemical  reactions.  The  experiments  indicated 

that  for  film  boiling  a  region  exists  for  which  a  stable 

vapor  layer  may  be  expected. 

Assuming 

-  steady,  two  dimensional  laminar  flow 

-  pressure  everywhere  constant 

-  heated  surface  temperature  constant 

-  temperature  of  liquid  vapor  interface  is  equal  to  satu¬ 
rated  liquid  temperature 

-  thermodynamic  and  transport  properties  of  liquid  and 
vapor  are  constant 

-  thermal  radiation  effect  on  vapor  and  liquid  layers  is 
neglected 

-  buoyancy  effect  on  forced  convection  boundary  layers  is 
neglected 

-  velocity  and  temperature  profiles  are  linear  in  vapor  and 
parabolic  in  liquid 


and  using  the  Pohlhausen  integral  method  the  results  of  the 
experiments  are  predicted.  In  some  of  the  experiments  dra¬ 
matic  instabilities  were  found,  and  the  authors  had  to  con¬ 
clude  that  producing  a  stable  vapor  layer  between  a  surface 
and  a  liquid  drastically  reduces  friction  drag  but  unstable, 
two-phase  flow  (like  nucleate  boiling  and  uncontrolled  sub¬ 
limation)  may  actually  cause  the  drag  to  increase. 


In  this  patent  Lang  "designed”  a  torpedo  capable  of 
using  the  gas  film  drag  reduction  phenomena.  The  lift  is 
provided  by  the  body  shape  and  the  exhaust  gases  are  used  to 
produce  the  gas  layer.  At  the  rear  end  of  the  torpedo  the 
gas  layer  is  removed  by  suction. 

10.  B.  R.  van  Driest  (111 

The  author  indicates  that  the  quest  for  higher 
speeds  of  underwater  craft  can  be  solved  by  using  the  knowl¬ 
edge  acquired  in  high  speed  aerodynamics. 

In  a  number  of  figures  the  author  gives  local  fric¬ 
tion  coefficients  as  a  function  of  Reynolds  number,  speed, 
wall  temperature  and  ocean  depth.  Further  effects  of  hydro- 
dynamic  heating,  cavitation  and  insertion  of  additives  are 
discussed. 

The  author  concludes  that  although  the  propulsion 
systems  have  to  be  developed,  it  is,  in  principle,  possible 
to  obtain  high  velocities  for  underwater  vehicles. 


B.  QUO  VADIS 

Considering  the  above  discussion  it  was  felt  that  more 
knowledge  had  to  be  obtained  of  the  phenomena  taking  place 
in  the  two-phase  boundary  layer.  A  first  step  would  be  to 
extend  the  work  on  flat  plates  to  two-dimensional  bodies 
with  pressure-gradients.  It  is  logical  to  follow  the  paper 
of  Sparrow,  Jonsson  and  Eckert  [7],  discussed  above  in  Sec¬ 


tion  A-6. 


III.  THEORY 


A.  GOVERNING  EQUATIONS 

To  predict  the  skin  friction  end  drag  char&cterics  from 
first  principles  the  starting  position  is  based  on  the  basic 
conservation  laws  governing  the  velocity  distribution,  here 
in  Cartesian  coordinates. 

Continuity;  |f  ♦  ♦  IttH  +  2j|!£L  .  0 


Eqn  of  State: 


+  &[v(ly  +  I*)]  +  A[y(lx  +  Is)] 

p5i  "  *  ‘  I?  +  w[|1(2ly  "  fdi'r  "]] 

*  &[u(l*  +  1y]]  +  ^[“(ly  +  In]] 

"m  ■  2  *  If +  &[u(2H  *  fdlv  5)] 

+  &[u(ls +  If)] +  w[u(lf  + !?]] 


0  Rg  * 


U  ■  W(PyT) 
dO  dW 

sssssl  aaation»  g  ■  *ar  +  at 

where  w  is  the  velocity  with  components  u,  v,  w  in  respec¬ 
tively  the  x,  y  and  z  direction,  is  the  operator 
3  3  3  3 

yfe  ♦  ♦  Vj—  +  ,  X,  Y  and  Z  are  the  body  forces  in 
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respectively  x,  y  and  z  directions,  P  is  the  local  pressure, 

Sg  is  the  gas  constant,  Q  is  the  heat  into  the  system,  ^  is 

the  internal  energy  and  W  is  the  work  done  by  the  system. 

Assuming  steady  state,  absence  of  body  forces,  viscosity 

3P 

y(T,P)  is  constant,  two-dimensional  flow,  ^  is  negligible, 
laminar  flow,  and  both  fluids  incompressible  (the  energy 
equation  decouples) ,  these  equations  simplify  to  the  Prandtl 
boundary  layer  equations: 

Continuity:  lx  +  ”ly  “  ® 

Navier-Stokes:  «|H  ♦  v|H  .  -  A  f£  + 

Using  the  subscript  g  for  properties  of  the  inner  fluid  and 
the  subscript  L  for  the  outer  fluid,  it  is  possible  to  write 
two  sets  of  equations  describing  the  two-fluid  laminar 
boundary  layer. 


Inner  fluid:  ■  0 


-3u  ,  „3u 


3P 

Tx 


32u 


(1) 

(2) 


Outer  fluid:  +  -|^ 


0 


(3) 


3u  .  3u  1  3P  .  32u 

*3*  +  *5y  m  '  tlx  +  VL^I 


(4) 


To  describe  the  problem  fully,  however ,  it  is  required  to 


state  the  boundary  conditions  and  the  conditions  at  the 
interface  between  the  two  fluids.  At  the  wall,  in  the  inner 


fluid,  the  velocity  parallel  to  the  wall,  u,  is  zero  (no 
slip)  and  the  velocity  perpendicular  to  the  wall  is  equal  to 
the  imposed  blowing  velocity.  At  the  interface  the  tangen¬ 
tial  velocity  and  shear  stress  have  to  be  continuous.  Fur¬ 
thermore,  the  mass  transfer  across  the  interface  is  assumed 
to  be  zero.  In  the  outer  fluid  u  has  to  approach  U(x)  as  y 
approaches  infinity.  A  formal  statement  of  these  conditions 
is: 

at  y  ■  0  u  »  0 

v  -  vw 

at  interface  u^  ■  u^ 

3u  3u_ 

Mijli  ‘  “iHi 

i  ■  0  or  v„  (interface)  *  v_  (interface)  ■  0 
g  I, 

in  outer  fluid  u  0(x)  as  y  » 

For  the  given  assumptions,  equations  1-4  and  the  boundary 
conditions  describe  the  problem  completely.  This  system  pre¬ 
sents  substantial  difficulties  in  attempting  a  general  solu¬ 
tion.  Similarity  solutions,  however,  will  solve  the  deadlock. 

B.  TRANSFORMED  EQUATIONS 

Similarity-type  solutions  are  well  known  and  have  the 
characteristic  that  the  velocity  profiles  have  similar  shapes 
at  all  x-positions  along  the  surface.  Mathematically  this 
gives  the  enormous  simplification  that  the  governing  partial 
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differential  equations  are  converted  into  ordinary  different 
tial  equations.  Before  defining  the  dimensionless  para¬ 
meters  it  should  be  noted  that  the  continuity  equations  are 
satisfied  by  the  stream-function  $  defined  by 

“  -  and  v  -  - 

Following  Schlichting  [12] ,  the  following  dimensionless 
parameters  can  be  introduced  for  the  inner  fluid: 

U  L 


Reynolds  number  Re 


(where  end  L  are  some 
reference  quantities) 


y/ssi  x 

ng  *  £g‘(xT  5g  *  £  (where  g(x)  is  a  dimensionless  scale 

factor  for  y) 

<Mx,y)  /I e~  r  u(x)  a(x) 

*<vv  •  t  rerSEf  or  *<x'*>  ■  flVV~ 

To  simplify  further  manipulation ,  it  is  wise  to  state  at 

this  point  that  to  ensure  similarity  f(n?»^)  can  not  depend 

on  -  E  (x)  so  that  f  *  f(n  ). 

9  9  9 

From  the  above,  the  velocity  components  may  be  derived, 
u  •  0(x)  f •  (rig) 

and  v  -  -  «<*>>  +  L°W-gtS>  ngf»'  'V  S,(x) 


z*®;  ”  z5®; 

Inserting  these  expressions  into  the  Prandtl  boundary  layer 
equation  and  after  a  considerable  amount  of  algebra  the  equa¬ 
tion  converts  into 
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I  I  I 


+  a  f  f  " 
<3  9 


+  8(1  -  f  •*)  «  0 


where  a  *  L  §--x)  gj(U(x)  g(x)) 
00 

6  .  *•  g2(*>  dD(x) 

"  0. 


Reasoning  in  the  same  way  for  the  outer  fluid  after 
defining 

•  «i-r 

Jj 

<Mx,y)  /5e£ 
fL<nL#?L)  *  L  uTxTgTxT 

one  arrives  at 


fL"'  +  afLfL"  +  8(1  -  fL'2)  -  0 


where  a  and  8  are  defined  the  same  as  for  the  inner  fluid. 

It  should  be  noted  that  to  satisfy  the  similarity  re¬ 
quirement  a  and  8  should  not  depend  on  x.  Still  following 
Schlichting  [12],  it  is  possible  to  arrive  at  the  following 
general  conditions  for  U(x)  and  g(x): 


EiiL 

U» 


where  K  is  an  integration  constant 


g(x) 


l2a  * 9)!  tni] 


Now,  without  loss  of  generality,  take  a  *  1  and  introduce  a 

8  U  (x) 

new  constant  m  -  so  that  the  expressions  for  and 

CO 

g(x)  convert  to 


I 


U(X) 

u_ 


g(x) 


(1+m) f  2  x 
i+m  L 


m 


2  x  U» 
m+T  L  UTxT 


Furthermore,  the  transformation  equations  for  the  ordinates 


are  ng  -  y 


/m+1  U  (x)“  , 

/t  "cnr  “d 


,x 


g-  ~  V  ~  L 

With  Schlichting  [12],  it  can  be  noted  that  these  re¬ 
quirements  can  be  satisfied  by  the  flow  in  the  neighbourhood 
of  the  stagnation  point  of  a  wedge  with  included  angle  Bit, 
where  U(x)  is  given  by  Cxm. 

Next  the  boundary  and  interface  conditions  have  to  be 
converted  by  expressing  them  into  the  newly  defined  dimen¬ 
sionless  parameters 
At  y  *  0 


or  n 


U(x)  fg* (n)  =  0  or  f  • (0) 


w 


-L  f  (n  )  . 

- 2—a-^{0(x)  g(x)> 


/rS 


or  fg<V 


"vw  /Reg 

,  d  (U  (x)  q  (x) ) 
L  - 


at  interface  n 


i,i  »  *  0  or  fg(ngi)  -  £l(0) 


\  or  V'V*  -  £l'(0) 


du 


dUL 


Ug  "W  “  UL  "3y  or  fg’ '  <  V  *  R  fL'  '  (0)  where  R  “ 
for  y  +  »  or  nT  -*•  ®  u  -*■  U  (x)  or  fT  '  -*■  1  as  nT  -*■  ®. 


(pu). 


(put 
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It  should  be  noted  that 

-  to  ensure  that  fL(n_)  is  independent  of  x: 

g  g 

vw  ^  )  .  Inserting  the  expressions  for  U(x) 

and  g(x)  one  obtains  vw  ^ 

-  R  depends  on  the  properties  of  both  fluids  depending  on 
temperature  and  pressure. 

Summarizing  the  problem  can  be  formally  stated  by 

equations:  f  ' ' '  +  f  f  ’ '  +  B(1  -  f'2)  *  0  (5a) 

g  g  g  g 

fL'"  +  fL  fL"  +  8(1  -  fg|2)  *  0 
and  conditions:  f  ' (0)  »  0 


(5b) 


IV.  EXPLORED  WAYS  OF  SOLUTION 


A.  SERIES  SOLUTION 

Following  the  procedure  of  Sparrow,  Jons son  and  Eckert 
[7],  a  series  solution  of  the  momentum  equation  seems  to  be 
appropriate.  Expanding  (n^)  with  respect  to  ■  0  gives 

’ 10) n„2  £' * 1 (0) n  3 

Vng>  ’  V0)  +  fg'(0,ng  +  -g"Ti  '  3  +  “a — 3! - *-  +  •  •  • 

The  first  two  terms  in  the  expansion  axe  available  from  the 
boundary  conditions.  The  f^‘ 1 (0)  is  unknown  for  now  and  will 
be  determined  by  the  solution.  Higher  derivatives  cam  be 
expressed  in  terms  of  lower  derivatives  by  taking  sequen¬ 
tially  derivatives  of  both  sides  of  the  momentum  equation 

£•..  ,.f  f.<  -  0(1  -  f'2) 

and  evaluating  at  rig  ■  0. 

Some  expressions  for  derivatives  are: 

f (III)  (0)  *  -f  (0)  f '  M0)  -  3 

f (IV)  (0)  »  f2(0)f '  '  (0)  +  8f  (0) 

f(V)(0)  -  -f3(0)f'M0)  -  Bf 2  (0)  +  (28-1)  f ' ' 2  (0) 

f (VI)  (0)  -  f4(0)f'M0)  +  Sf 3  (0)  -  (83-5)  f  (0)  f  *  * 2  (0) 

-  (68-4)  Sf '  M0) 


It  was  impossible  to  find  a  recurrency  in  the  expressions. 

On  applying  the  interface  condition  f ( n  *  0  makes  it  pos¬ 

sible  to  write 


It  is  now  possible  to  express  f  ' 1  (0)  in  liquid  and  gas  para¬ 
meters.  Here  Sparrow,  Jonsson  and  Eckert  refer  to  an  earlier 
paper  by  Cess  and  Sparrow  [3] .  This  paper  is  discussed  in 
Chapter  II,  Section  A3.  This  paper  gives  a  table  relating 
f^' (0)  and  fL' ' (0)  which  corresponds  to  a  solution  of  the 
equation  for  a  flat  plat. 

At  this  point  it  was  realized  that  the  data  necessary  to 
go  any  further  was  not  available,  therefore  this  method  was 
abandoned,  and  other  methods  of  solution  were  explored. 

B.  NUMERICAL  SOLUTIONS 
1.  IODE  [13] 

IODE  is  an  interactive  ordinary  differential  equation 
solver  developed  by  the  Computer  Center  of  the  Naval  Post¬ 
graduate  School.  The  program  accepts  a  number  of  first  order 
non-linear  ordinary  differential  equations  with  associated 
initial  conditions  and  solves  the  equations. 

In  this  case  equations  5  were  converted  into : 


dy2 

“2F  “  y3  "  f"  (6b) 

dy3  2 

m  £•••  m  -y1y2  -  3  +  sy22  (6c) 

where  y^,  y2  and  y3  are  dummy  variables  for  the  derivatives 
of  f(n)#  with  initial  conditions 

yx (0)  -  f(0)  (7a) 

y2  (0)  *  f  (0)  (7b) 

y3(0)  -  f"(0)  (7c) 

The  equations  6  are  the  same  for  both  fluids  but  as  noted 

before  the  initial  conditions  differ. 

Sanders  [14]  suggested  that  it  should  be  possible  to 

guess  the  thickness  of  the  inner  fluid  layer  ngi  and  then  a 

value  of  fg,( (0)  could  be  found  by  trial  and  error  so  that 

no  mass  transfer  occurs  across  the  interface  (f  (n  j)  ■  0) . 

g  gi 

The  values  for  f  .  f '  and  f'*  at  the  interface  give  the 

g  9  9 

initial  conditions  for  the  outer  fluid  equations  6.  The 
values  of  f^' ( n^)  should  approach  unity  for  large  enough  nL 
while  fL''(nL)  approaches  0.  If  ^'(n^)  does  not  go  to  unity 
the  initial  guess  for  the  layer  thickness  was  wrong  and 
another  value  for  the  inner  fluid  layer  thickness 
chosen.  This  process  is  repeated  until  a  layer  thickness  is 
found  for  which  f^'  («)  ■  1  and  f^'M*)  “0. 

The  calculations  carried  out  on  the  IBM  3033AP  com¬ 
puter  system  showed  that  this  procedure  worked ,  but  while  the 
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calculations  tiara  fast  tha  intaractiva  sida  of  tha  program 
has  provan  to  ba  too  tima  consuming. 


2.  Dash  Top  Computer 

So  tha  necessity  arose  to  find  soma  way  to  decrease 
tha  tima  spent  communicating  with  tha  computer  and  to  auto¬ 
mata  the  trial  and  error  process.  Furthermore,  it  was 
necessary  to  find  sons  way  to  interrupt  tha  execution  and 
to  execute  the  program  by  steps  to  determine  instabilities 
in  the  solution. 

Since  a  Wang  2200  desk-top  computer  with  line  printer 
was  available,  the  problem  was  rewritten  for  this  computer 
and  the  solution  for  a  number  of  cases  was  carried  out.  It 
has  to  be  realized  that  a  desk-top  computer  is  relatively 
slow  so  that  a  compromise  between  calculation  time  and  ac¬ 
curacy  of  the  result  had  to  be  found.  It  was  decided  that 
the  primary  goal  was  to  find  *a  method  of  solution  and  not  to 
tabulate  a  large  amount  of  accurate  data.  In  the  light  of 
all  the  assumptions  high  accuracy  would  be  doubtful  anyway. 


V.  SOLUTIONS 


A.  DEFINITION  OF  THE  PROBLEM 

The  momentum  equation  will  be  written  as  three  first 
order  differential  equations  as  for  the  MIODE  program": 

dyl 

-ar  -  *2  •  f‘ 

dy2 

-a#  ■  *3  * 

dy3  2 

-ar  -  -8-^3 

Now  multiplying  both  sides  of  these  equations  by  dn  and 
changing  to  finite  differences  one  arrives  at 

Ayx  -  Y2*An 

Ay2  "  y3*An 

ay3  *  <-yx*y3  -  s  -  sy12)*An 

Now  by  stepping  n  through  the  desired  range  of  values, 
starting  with  known  or  guessed  initial  values  for  n^,  ,  y2 

and  y3  it  is  possible  to  calculate  y^,  y2  and  y3  as  a  func¬ 
tion  of  n. 

It  has  to  be  noted  at  this  point  that  we  approximate  the 
integrals  by  a  "Lower  Riemann  Sum".  It  is  known,  however, 
that  using  a  "sufficient"  number  of  steps  guarantees  a  satis¬ 
factory  accuracy.  This  procedure  is  applicable  for  both  the 
inner  and  the  outer  fluid.  The  initial  values  for  the  outer 


fluid  follow  logically  from  the  solution  of  the  inner  fluid 
and  the  interface  conditions. 

Given  values  for  the  blowing  parameter  f^(0)  and  the 
properties  of  the  fluids  parameter  R  the  procedure  is  as 
follows: 

-  Choose  an  inner  fluid  layer  thickness. 

-  Find  f''(0)  by  an  interval  halving  method  so  that  f(n  ,) 
(no  mass  transfer  across  the  interface)  equals  0,  or  91 
in  our  program  is  smaller  than  a  chosen  value. 

Note:  f  (0)  ■  given;  f '  (0)  ■  0  and  f  "  (0)  is  chosen  as 
indicated. 

-  Apply  the  interface  conditions  on  the  values  for  f,  f' 
and  f'1  at  the  interface  to  find  the  initial  values  for 
the  solution  of  the  outer  fluid  equations.  Choose  an 
end  value  for  nL  large  enough  to  be  well  "outside”  the 
boundary  layer. 

-  Check  if  fr  '  (*)  *  1  but  also  if  f  ' '  (°°)  a  0  as  it  is 
possible  that  an  unstable  solution  is  found.  In  case 
the  conditions  are  satisfied  the  solution  is  found.  If 
not  another  guess  on  the  inner  fluid  layer  thickness  has 
to  be  made. 

-  The  value  of  the  inner  fluid  layer  thickness  is  chosen 
in  another  interval  halving  process. 

•  This  process  is  repeated  until  the  requirements  of 
f  ’ (®)  ■  1  and  fL' ' (•)  a  0  are  met. 

A  copy  of  this  program  is  enclosed  as  Appendix  A. 

From  the  above ,  it  follows  that  the  initial  intervals 

for  f  'MO)  and  ri  <  have  to  be  determined.  The  interval  for 
g  gi 

ngi  should  be  as  narrow  as  possible  and  the  interval  for 
f''(0)  should  be  a  little  wider  than  expected  to  prevent  in¬ 
stabilities  and  "limit-riding".  These  values  were  determined 
by  trial  and  error  and  the  results  are  discussed  in  a  later 
chapter. 


B.  TESTING  THE  PROGRAM 

1.  Number  of  Steps  of  Independent  Variable 

The  program  is  used  with  100  steps  of  the  independent 
variable.  A  sample  problem  was  solved  with  both  100  and  200 
steps.  The  deviations  were  in  the  fourth  decimal  (see  Table 
5.1).  The  time  required  for  the  calculation  however  was 
doubled  from  approximately  45  minutes  to  90  minutes. 

2.  One  Fluid  Flow 

It  is  clear  that  when  R  approaches  unity,  the  two 
fluids  are  hydrodynamically  identical.  Computer  solutions 
were  obtained  for  a  single  fluid  flow  for  0*0  (flat  plate) 
and  0*1  (stagnation  flow)  and  further  for  0  *  0.2,  0.4, 

0.6  and  0.8.  The  results  were  in  good  agreement  with  Nickel 
[5],  "abbildung  8".  (See  Table  5.2.) 

3.  Solution  for  R  Close  to  1 

A  solution  with  a  R  of  1.1  was  compared  to  that  for 
a  one-fluid  problem  and  the  results  were  close  enough  to 
give  further  confidence  in  the  correctness  of  the  program  to 
handle  non-unity  values  of  R.  (See  Table  5.3.) 


4.  Solution  for  Small  Wedge  Angles 


VZ.  RESULTS 


A.  DEFINITION  OF  A  SAMPLE  PROBLEM 

To  demonstrate  tha  capabilities  of  tha  program,  a  spacial 
caaa  was  chosen  for  thorough  investigation.  For  this  case, 
"air  at  standard  conditions"  and  watar  war a  chosan  as  tha 
two  fluids  so  that  R  ■  200.  For  a  range  of  blowing  para¬ 
meters  fg(0),  and  wedge  angles  8,  tha  velocity  profile  was 
calculated  and  tha  layer  thickness  ngi  of  tha  inner  fluid 
and  the  shear  stress  parameter  fg* ' (0)  at  the  wall  were 
determined. 

B.  THICKNESS  OF  THE  INNER  FLUID  LAYER 

Figure  6.1  gives  the  dimensionless  thickness  of  the  inner 
fluid  layer  ngi  as  a  function  of  the  blowing  parameter  fg(0) 
for  various  wedge  angles  8.  Figure  6.2  is  an  enlargement  of 
these  curves  for  small  blowing-parameters.  These  graphs  lead 
to  the  following  observations s 

-  por  small  values  of  the  blowing  parameter  the  curves  are 
straight  lines,  Hgi  decreases  very  little  with  increasing 
wedge  angle  and  if  directly  proportional  to  fg(0). 

-  For  larger  blowing  parameters  fg(0),  r\„^  increases  more 
rapidly  with  increasing  fg(0).  For  lafge  blowing  para¬ 
meters  and  small  angles  the  thickness  of  the  inner  fluid 
layer  increases  drastically  with  increased  blowing 
parameter . 

-  For  large  values  of  the  blowing  parameter,  the  thickness 
of  the  inner  fluid  layer  decreases  rapidly  with  increas¬ 
ing  wedge  angle  8. 
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C.  WALL  SHEAR  STRESS 

Figures  6.3  and  6.4  give  the  dimensionless  wall  shear 

stress  f  ' '  (0 )  as  a  function  of  the  blowing  parameter  f(0) 

9 

and  the  wedge  angle  8. 

-  For  all  values  of  wedge  angle  8#  blowing  decreases  the 
wall  shear  stress  fg'MO).  This  reduction  is  spectacular 
for  small  values  of  the  blowing  parameter;  e.g.  for  both 
the  flat  plate  (6  ■  0)  and  stagnation  flow  (8*1) 

f-' * (0)  is  reduced  to  half  the  non-blowing  value  when 
fg(0)  -  0.0025. 

-  For  all  values  of  the  blowing  parameter  the  wall  shear 
stress  f'(0)  increases  with  increasing  wedge  angle  8; 
the  smaller  thickness  of  the  inner  fluid  layer  caused  by 
the  higher  potential  velocity. 

-  For  large  values  of  the  blowing  parameter,  the  wall  shear 
stress  is  very  small  and  nearly  equal  for  all  wedge 
angles.  The  values  obtained  indicate  that  the  wall  shear 
stress  for  large  blowing  parameters  approaches  the  wall 
shear  stress  for  a  wedge  in  an  infinite  volume  of  inner 
fluid.  (See  Table  6.1.) 


D.  VELOCITY  PROFILES 

Figure  6.5  gives  velocity  profiles  (f ' ( n)  vs.  n)  for 

various  blowing  parameters  for  a  wedge  angle  8  *  0.2  (36°). 

The  and  have  been  expressed  in  the  same  relative  scale 

y/Ke^  (y-<5)/Re~ 

by  noting  that  *  L  and  nL  ■  1  L’  ~g(X) —  resulting  in 


The  following  observations  can  be  made: 

-  For  the  smaller  blowing  parameters  the  velocity  profiles 
in  the  inner  fluid  are  nearly  straight  lines. 

-  The  stronger  the  blowing  the  thicker  the  layer  of  the 
inner  fluid  causing  a  larger  proportion  of  the  potential 
velocity  to  be  obtained  in  the  inner  fluid  layer.  For 
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small  blowing  parameters  this  proportion  of  the  potential 
velocity  increases  rapidly  with  increasing  blowing  para¬ 
meter.  Figure  6.6  gives  the  proportion  of  the  potential 
velocity  (or  f  ' (nai))  as  a  function  of  the  blowing  para¬ 
meter  fg(0)*  9  9 

-  Particularly  interesting  is  the  fact  that  the  velocity 
profiles  in  the  inner  fluid  and  the  outer  fluid  have  op¬ 
posite  curvatures.  This  is  comparable  to  the  boundary 
layers  associated  with  the  gas-filled  cavity  attached  to 
a  body  after  water  entry.  The  cavity  around  such  a  body 
has  a  stable  interface;  breakdown  occurs  only  when  too 
much  gas  is  left  at  the  rear  separation.  Although  no 
rigorous  proof  can  be  given  it  is  felt  that  this  charac¬ 
teristic  of  the  velocity  profile  provides  stability  of 
the  interface  in  two  fluid  boundary  layer  flow. 


E.  DRAG  REDUCTION 

To  characterise  the  amount  of  drag  reduction,  the  local 


drag  parameter  f  ' ' (0)  obtained  for  the  two  fluid  flow  was 


compared  to  the  local  drag  parameter  obtained  without 
blowing. 

To  compare  these  local  drags,  a  minor  converson  has  to 
be  made.  The  shear  stress  at  the  wall  is 


■*(£) 


y»0 


In  terms  of  the  transformed  variables,  the  expression  for  t 
becomes 

U  (x) 


*  -  vg  EgTxf  /®*g  <*"<">  >n-0 
which  can  be  expressed  in  the  outer  fluid  parameters 


For  the  outer  fluid,  without  blowing,  the  reference  shear 


stress  Tq  is: 

*0  <V  <">>n-0 

fg(0)-0 


Now  taking  the  ratio: 


T 


fa"(0) 

^ri^tgi0) 


0 


The  values  of  this  reference  shear  stress  parameter  f  ' '(0) 
are  calculated  and  reflected  in  Figure  6.7.  Since  R  is  a 
known  constant,  for  each  value  of  the  shear  stress  parameter 
fg''(0)  it  is  possible  to  determine  the  ratio  t/Tq.  Table 
6.2  gives  the  values  of  t/tq  as  a  function  of  wedge  angle  S 
and  blowing  parameter  fg(0)  for  R  *  200. 

Inspection  of  this  table  reveals  that  a  drag  reduction 
is  observed  for  all  blowing  parameters  and  wedge  angles.  It 
is  remarkable  that  for  small  values  of  the  blowing  parameter 
the  ratio  t/Tq  decreases  for  increasing  wedge  angle  8  while 
for  large  values  of  the  blowing  parameter  f(0)  the  ratio 
t/Tq  increases  with  increasing  wedge  angle  S.  For  all  wedge 
angles  3,  however,  the  ratio  t/tq  decreases  with  increasing 
blowing  strength. 


F.  REQUIRED  POWER 

The  power  required  to  maintain  blowing  is 

Power  ■  P  •  vw 

where  P  is  the  local  pressure  given  by 


P  “  PSTAG  +  3x  ^ 


where  PSTAG  i*  the  stagnation  pressure  at  the  leading  edge 
of  the  wedge  and  ||  -  -pgUgj[  »-Pg  C2  m  x2®*1  <0.  If  the 
inner  fluid  is  at  stagnation  pressure  PgTAG  a  pressure  dif- 
ference  of  at  the  wall  exists  so  that  the  fluid  flows  out 
by  itself.  The  flow  rate  then  can  be  controlled  by  the 
blowing  hole's  size.  Of  course,  the  ejected  fluid  has  to  be 
replaced  requiring  it  to  be  brought  up  to  PSTAG*  In  practice 
this  is  of  less  interest  as  a  considerable  amount  of  exhaust 
gases  of  the  propulsion  systems  at  high  pressures  should  be 
available . 


II.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  has  been  shown  theoretically  that  a  dramatic  drag 
reduction  is  possible  for  the  infinite  wedge  by  blowing  with 
a  lighter  and  less  viscous  fluid,  and  that  the  parameters 
involved  can  be  calculated  with  relatively  simple  means. 
Furthermore  the  developed  method  gives  another  way  of  solving 
the  flat  plate  case  and  the  stagnation  flow  case. 

This  work  directly  invites  a  number  of  follow-on  investi¬ 
gations  of  phenomena  associated  with: 

-  finite  two-dimensional  bodies 

-  three-dimensional  bodies 


-  a  compressible  inner  fluid 

-  experimental  verification  of  the  theory 

-  stability  of  the  interface 
all  of  equal  importance. 


APPENDIX  A 


"BASIC"  PROGRAM  FOR  WANG  DESK  TOP  COMPUTER 
TO  SOLVE  TWO  FLUID  LAMINAR  FLOW  CHARACTERISTICS 


1  REM  THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  F,F',F" 

2  REM  AND  NCI  FOR  A  TWO-LIQUID-LAMINAR-BOUNDARY-LAYER 

3  REM  GIVEN  A  BLOWING -PARAHETER  F<0) , A  WEDGEANGLE  B  AND 

4  REM  A  LIQUID-PROP ERTIES-CONSTANT  R. 

5  REM 

6  REM  TO  USE  THE  PROGRAM  ONE  HAS  TO  GUESS  ON  AN  INITIAL 

7  REM  INTERVAL  FOR  THE  I NNER-FLU ID-LAYER -THICKNESS  AND 

8  REM  FOR  THE  LOCAL  SHEARS7RESS-PARAMETER  F" (O). THESE 

9  REM  VALUES  CAN  BE  OBTAINED  FROM  GRAPHS  IN  THE  THESIS. 

10  REM  ONE  HAS  TO  TAKE  THE  INTERVAL  FOR  THE  THICKNESS  AS 

11  REM  SMALL  AS  POSSIBLE  AND  THE  INTERVAL  FOR  THE  SHEAR- 

12  REM  STRESS  A  LITTLE  BIT  LARGER  THAN  EXPECTED  TO  PREVENT 

13  REM  INSTABILITIES  AND  LIMITRIDING. 

14  REM 

15  REM  IN  SHORT  INSERT  VALUES  FOR 

16  REM  -  F(0)  AS  Cl <1>, LOWER  LIMIT  OF  F" <0>  AS  G6  AM) 

17  REM  -  UPPER  LIMIT  OF  F*  <0>  AS  G7  IN  LIKE  120 

18  REM  -  LOWER  LIMIT  OF  NCI  AS  H6  AND  UPPER  LIMIT  OF 

19  REM  -  NGI  AS  H7  IN  LINE  130 

20  REM  -  WEDGE  ANGLE  B  AND  FLUID  PROPERTIES  CONSTANT  R 

21  REM  -  IN  LINE  140 

22  REM 

23  REM  THE  OUTPUT  WILL  CONSIST  OF  TWO  TABLES  GIVING 

24  REM  FG,FG ' ,FC" ,FL, FL'  AND  FL"  AS  A  FUNCTION  OF  THE 

25  REM  DIHENSIQM.ES  PARAMETERS  NG  AM)  NL.PLUS  A  HEADING 

26  REM  OF  THE  MOST  IMPORTANT  PARAMETERS  AND  CONSTANTS. 

27  REM 

28  REM 

100  DIM  G0(101),  Gl(lOl),  G2< 101 ) ,  G3(101> 

105  DIM  L0< 101 ) .  LI (101),  L2<101),  L3(101) 

110  G0<1)*0:  G1 < 1 )*-. 5:  G2(l)«0:  G6-0:  G7»2 

115  H6-.95:  H7-1.05 

120  G4-G6:  G5*G7 

125  A>1:  8-1:  R-200 

130  Z-0:  Zl-0 

135  H8-H6:  H9-H7 

140  H»<H8H-f9)/2 

145  GS-G6:  G9-G7 

ISO  C3<l)*<G8+G9)/2 

155  IF  H-0  THEN  400 

160  IF  Z  <  35  THEN  180 


16S  PRINT  “AFTER  35  ITERATIONS  IN  GAS  UNABLE  TO  >t2T  REG'S 
170  PRINT  H, G1 < 101 ) , G2< 101),G3<1) 

175  GOTO  495 
180  HO-H/lOO 
185  FOR  1-1  TO  100 
130  H1-C2(I)*H0 

195  H2«G3(I)«H0 

200  H3*H0#(B#G2< I )*C2< I )-B-A*Ci ( I >*G3( I ) ) 

205  G1(1*1)>G1<I)+H1 

210  G2(I+1>*G2<I)*H2 

215  G3<I*1)»G3<I)+H3 

220  G0(1«-1>-C0(I>+H0 

22S  NEXT  I 
230  Z-Z+l 

235  IF  ABS<C1 ( 101 > )  <  .00000001  THEN  255 
240  IF  Cl (101)  >0  THEN  250 
24S  G8»C3(1>:  GOTO  ISO 

250  C9«G3<1>:  GOTO  150 

255  L0(1)«0:  L1(1)-G1(101):L2(1>-G2(101>:  L3( 1 )»C3< 101 )/R 

250  N»10 

255  M0-M/100 

270  FOR  1-1  TO  100 

275  M1-L2(I)*M0 

280  M2»L3<I)*M0 

285  M3«MO*<B*L2(I)*L2(I>-B-A«Ll(I>*L3(I>> 

290  L1<X+1>"L1<I)+M1 

295  L2<  I+l  )*L2(  I )  ♦M2 

300  L3<X+1>»L3(I>*M3 

305  L0< I+l )»L0< I )>M0 

310  NEXT  I 
315  Zl-Zl+1 

320  IF  Zl<  35  THEN  335 

325  PRINT  “AFTER  35  ITERATIONS  IN  LIQ  UNABLE  TO  MEET  REG'S 

330  GOTO  415 

335  C5»L2(101)-1 

340  IF  H-O  THEN  370 

345  IF  ABS<L3( 101 ) )  >  .0001  THEN  355 

350  IF  ASS<C5>  <  .00001  THEN  410 

355  IF  CS  >  O  THEN  355 

350  H9*H:  G6»G3(1):  Z»0:  GOTO  140 

365  H8»H:  G7-C3(l>:  Z-0:  GOTO  140 

370  IF  ABS<L3( 101 ) )  >  .0001  THEN  380 

375  IF  ABS(C5>  <  .00001  THEN  410 

380  IF  CS  >  0  THEN  390 

385  G8*G3(1):  GOTO  395 

390  G9«G3(1>:  GOTO  395 

395  G3<l)>(G8+G9)/2 

400  L3( 1 )»G3< 1 ):  L0( 1 )*0:  L1<1)«G1<1>:  L2(1)*L2(1) 

405  GOTO  250 

410  PRINT  "LAMINAR  FLOW  OVER  WEDGE  WITH  BLOWING" 

415  PRINT  :  PRINT  :  PRINT 


.*.  .  s  .%  /.  /-V 


PRINT  *8  ■  ■ ,B 

PRINT  *R  *  •  p 

PRINT  -FG(0>  =  %G1<1),“ INITIAL  LIMITS 

PRINT  "FG0PRIME(0>*  - , G3( 1 ) , G4,G5 
PRINT  *NI  »  * , G0< 101 ) , H6,H7 

PRINT  :  PRINT  :  PRINT 

PRINT  “NG" , "FC" , "FGPRIME* , ■FGDPRIME" :  PRINT 
FOR  1*1  TO  101  STEP  10 

PRINT  G0< I ) , G1 < I ) , G2< I ) , G3< I ) 

NEXT  I 

PRINT  !  PRINT  !  PRINT 

PRINT  "NL" , "FL* , "FLPRIME* , *FLDPRIME" :  PRINT 
FOR  1*1  TO  101  STEP  10 

PRINT  L0< I ) , LI < I ) , L2< I ) , L3< I ) 

NEXT  I 

PRINT  :  PRINT  :  PRINT  :  PRINT  "  ■  -  ■  •“£»*>• 
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Figure  6.3  Dimensionless  wall  shear  stress  as  a 

function  of  blowing  parameter  and  wedge 


Dimensionless  wall  shear  stress  for  the 
outer  fluid  without  blowing  as  a  function 
of  the  wedge  angle 
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0 

m 
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0.0585 
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0.0551 

0.0476 

0.0426 

0.0389 
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0.0468 
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0.0325 
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0.0266 

0.1 

0.0526 

0.0355 
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0.15 

0.0353 
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0.0167 
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0.0138 
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0.35 
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0.00692 

0.00656 
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0.00502 

0.75 

0.00337 

0.00351 

0.00363 
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0.00379 

0.00386 

1.0 
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0.00212 

0.00261 

0.00289 

0.00309 

0.00324 
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